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ABSTRACT 



SOLERI 6 , R. WlJNANDS 1 



We report on the discovery and the timing analysis of the first eclipsing accretion-powered mil- 
lisecond X-ray pulsar (AMXP): SWIFT J1749. 4-2807. The neutron star rotates at a frequency of 
~517.9 Hz and is in a binary system with an orbital period of 8.8 hrs and a projected semi-major 
axis of ~1.90 lt-s. Assuming a neutron star between 0.8 and 2.2 Mq and using the mass function of 
the system and the eclipse half-angle, we constrain the mass of the companion and the inclination 
of the system to be in the ~0. 46-0. 81 M and ~ 74.4° — 77.3° range, respectively. To date, this 
is the tightest constraint on the orbital inclination of any AMXP. As in other AMXPs, the pulse 
profile shows harmonic content up to the 3rd overtone. However, this is the first AMXP to show a 1st 
overtone with rms amplitudes between ~ 6% and ~ 23%, which is the strongest ever seen, and which 
can be more than two times stronger than the fundamental. The fact that SWIFT J1749. 4-2807 is 
an eclipsing system which shows uncommonly strong harmonic content suggests that it might be the 
best source to date to set constraints on neutron star properties including compactness and geometry. 
Subject headings: pulsars: general — pulsars: individual (WIFT J1749. 4-2807) — stars: neutron 
binaries: eclipsing 



1. INTRODUCTION 

The first accreting millisecond X-ray pulsar (hereafter 
AMXP) was discovered in 1998 (SAX J1808.4-3658, see 
IWijnands fe van der Klisl I1998T ) and since then, a to- 
tal o f 13 AMXPs h ave been found and studied in de- 
tail (|Patrunol I2010D . Most AMXPs show near sinu- 
soidal profiles during most of their outbursts. This 
is consistent with a picture in which only one of the 
hotspots (at the magnetic poles) is visible (see ref. be- 
low). Deviations from a sinusoidal profile (i.e, an in- 
crease in harmonic content) are generally interpreted as 
being caused by the antipodal spot becoming visible, per- 
hap s as accretion rate falls and the disk retreats (see, 
e.g. IPoutanen fe Gi crliriski 2003; Ibr agimov fc Poutanenl 
I2009L and references therein). 

Although the amplitude of the 1st overtone may reach 
that of the fundamental l ate in the outburst (see, e.g 
lHartman et~alll2008l[2009] ). no AMXP so far has shown 
pulse profiles where the 1st overtone is generally stronger 
than the fundamental throughout the outburst. 

The stability of the pulse profiles in some of the 
AMXPs means that pulse profile modeling can be used 
to set bounds on the compactness of the neutron star 
and hence the dense matter equation of state (see, e.g. 
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iPoutanen fc Gierlihskil [20031 : iPoutanen et alJl2009L and 
references therein). Unfortunately, there is often a large 
degeneracy between the parameters due to the number 
of free parameters needed to construct the model pro- 
file. One of these parameters is the inclination of the 
system, which to date has not been well-constrained for 
any AMXP. 

In this Letter we report on the discovery and tim- 
ing of the accretion-powered millisecond X-ray pulsar 
SWIFT J1749.4-2807. Thanks to the observed eclipses 
(jMarkwardt et all [2010). we set the tightest constraint 
on system inclination for any AMXP. This, coupled 
with the fact that the amplitude of the first overtone is 
higher/comparable to that of the fundamental for much 
of the outburst and that the amplitude of the first over- 
tone is unusually high, allows to put tight constraints on 
pulse profile models. We show that SWIFT J1749.4-2807 
has the potential to be one of the best sources for this 
approach to constraining the neutron star mass-radius 
relation and hence the EoS of dense matter. 

2. SWIFT J1749. 4-2807 

SWIFT J1749.4-2 807 was discovered in June 2, 2006 
(|Schadv et al.ll200"6T ). when a bright bur st was detected 
by the Swift burst alert telescope (BAT). IWijnands et al.l 
(2009) presented a detailed analysis of the Swift/BAT 
and Swift/XRT data and showed that the spectrum 
of the 2006 burst was consistent with that of a ther- 
monuclear Type I X -ray burst (jPalmer et al.l 120061 : 
iBeardmore et al.1l2006l see, also) from a source at a dis- 
tance of 6.7 ± 1.3 kpc. 

SWIFT J1749.4-2807 was detected again between 
April 10th an d 13th, 2010 using INTEGRAL and Swif t 
observations (|Pavan et al.l l2010t IChenevez et al.l 2010). 
We promptly triggered approved RXTE observations on 
this source to study X-ray bursts and to search for 
millisecond pulsations (Proposal ID:93085-09, PI: Wij- 
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Fig. 1. — Top panel: 2-10 keV flux as measured from 
RXTE/PCA and Swift/XRT observations. The flux of the last 
PCA observation (MJD 55307.5) is not shown; the spectrum of 
this observation was used as an estimate background emission (see 
text). Upper limits are quoted at 95% confidence level. We cal- 
culated the flux during the last PCA and last Swift/XRT detec- 
tion using WebPIMMS (assuming a power law spectrum with index 
1.8). Middle panels: Fractional rms amplitude and 95% confidence 
level upper limits of the fundamental and three overtones as a func- 
tion of time. Detections (> 3<j single trial) and upper limits are 
from ~ 500 and ~ 3000 sec datasets, respectively. Bottom panel: 
Ratio between the fractional rms amplitude of the 1st overtone and 
fundamental. Blue triangles represent points in which both har- 
monics are significantly detected in fti 500 second datasets. Black 
circles represent the ratio between the fractional rms amplitude of 
the 1st overtone and the 95% confidence level upper limit to the 
amplitude of the fundamental. Grey circles represent the same ra- 
tio but when the fundamental is significantly detected and not the 
1st overtone. This means that black circles represent lower limits 
while grey circles are upper limits. These ratios are independent 
of the background. 



nands). The first RXTE observation was performed on 
April 14th and lasted for about 1.6 ksec. We found strong 
coherent pulsations at « 517.9 Hz and at its first over- 
tone « 1035.8 Hz, showing that SWIF T J1749.4-2807 is 
an ac creting millisecond X-ray pulsar (|Altamirano et ahl 
l2010h . RXTE followed up the decay of the outburst on 
a daily basis. Preliminary results on the rms amplitude 
of the pulsations, orbital solution, discovery of eclipses, 
evolution of the outburst and upper limits on the quies- 
cent luminosity were rep orted in Astronom e r's Telegrams 
(lAltamirano et alj|2010l iBozzo et al.ll2010l; iBelloni et al. 



2010; Strohmayer & Markwardt 2010; Markwardt et al 



2010HYang et al.ll2010H Chakrab arty et al.ll2010f >. No op 



tical counterpart has been ident ified as yet, with a So- 
lower limit in the i-band of 19.6 (|Yang et al.l l2010). 

3. OBSERVATIONS, SPECTRAL ANALYSIS AND 
BACKGROUND ESTIMATION 

We used data from the Rossi X-ray Timing Explorer 
(RXTE) Proportion al Counter Array (P CA, for instru- 
ment information see lJahoda et al.ll2006| ). Between April 
14th and April 21st there were 15 pointed observations of 
SWIFT J1749.4-2807, each covering 1 or 2 consecutive 
90-min satellite orbits. 

We also analyzed data f rom Swift's X-ray telescope 
(XRT: iBurrows et al.l l2005h . There were a total of 10 
observations (target ID 31686), all obtained in the Pho- 
ton Counting (PC) mode. 

We used standard tools and procedures to extract en- 
ergy spectra from PCA Standard 2 data. We calculated 
response matrices and ancillary files for each observa- 
tion using the FT00LS routine PCARSP V10.1. Back- 
ground spectra were estimated using the faint-model in 
PCABACKEST (version 6.0). For the XRT, we used stan- 
dard procedures to process and analyze the PC mode 
data 7 . When necessary, an annular extraction region 
was used to correct for pile-up effects. We generated 
exposure maps with the task XRTEXP0MAP and ancillary 
response files were created with XRTMKARF. The latest re- 
sponse matrix files (v. 11) were obtained from the CALDB 
database. 

We used an absorbed power-law to fit all PCA/XRT 
observations. We first fitted all XRT spectra and found 
an average interstellar absorption of 3.5 x 10 22 cm~ 2 ; 
Njj was fixed to this value when fitting all (standard) 
background-subtracted PCA spectra. When comparing 
the fluxes estimated by PCA and XRT we found that 
the PCA fluxes were systematically higher. Only on 
one occasion RXTE and Swift observations were per- 
formed simultaneously (MJD 55306.69, i.e. at the end 
of the outburst) and in this case the flux difference was 
w 1.4 x 10 -10 erg cm -2 s _1 . This is consistent with that 
seen a day later. Since (i) the count rates during the last 
RXTE observation (when the source was below the PCA 
detection limit but detected by Swift/XRT) are consis- 
tent with those we measure during the eclipses (see Sec- 
tion [6|) and (ii) these count rates are consistent with the 
offset we find between PCA and XRT, we conclude that 
there is an additional source of background flux in our 
PCA observations. To correct for this, we also use the 
background-corrected spectrum of the last PCA obser- 
vation (OBSID: 95085-09-02-08, MJD 55307.5) as an es- 

7 http://www.swift.ac.uk/XRT.shtml 



TABLE 1 
Timing parameters for the AMXP SWIFT J1749. 4-2807 

Parameter Value 

Orbital period, P orb (days) 0.3673696(2) 

Projected semi major axis, a^sinz (lt-s) 1.89953(2) 

Time of ascending node, T asc (MJD) . . 55300.6522536(7) 

Eccentricity, c (95% c.l.) < 5.2 x 10~ 5 

Spin frequency v (Hz) 517.92001395(1) 

Pulsar mass function, f x (M ) 0.0545278(13) 

Minimum companion mass , M c (Mq) . 0.5898 

All errors are at Ax — 1. 
: The companion mass is estimated assuming a neutron star of 1.4 

M . 

timate of the additional source of background flux. This 
approach is optimal in crowded fields near the Galactic 
plane, where the contribution from the Galactic ridge 
emission and other X-ray so urces in the 1° PC A FoV 
becomes important (see, e.g.. lLinares et al.ll2007l l2008h . 

3.1. Background estimates and the fractional rms 
amplitudes 

Given the low flux during our observations, it is very 
important to accurately estimate the background emis- 
sion before calculating the pulse fractional rms ampli- 
tudes. The fact that the extra source of background 
photons is unknown complicates the estimation of total 
background flux as a function of time. For example, the 
background flux could be intrinsically varying; even in 
the case of a constant distribution of background flux in 
the sky, it would be possible to measure flux variations if 
the collimator (i.e. PCA) orientation on the sky changes 
between observations. Given the background uncertain- 
ties, we arbitrarily adopt as total 3-16 keV background 
per observation the modeled background plus a constant 
offset of « 17.5 ± 2 ct/s/PCU. This takes into account 
the w 19.5 ct/s/PCU as estimated by the eclipses, the 
last PCA observation and the PCA-XRT offset (which 
is equivalent to « 18 - 19 ct/s/PCU in the 3-16 keV 
band as estimated with WebPIMMS 8 and the best fit 
model to the XRT data), and the f=a 15.5 offset we would 
obtain if the additional source of background photons 
could change by ~ 20%. This conservative adopted pos- 
sible background range results in conservative errors on 
the pulsed fractions we report, i.e. the errors are proba- 
bly overestimated. 

4. OUTBURST EVOLUTION 

In Figure [T]we show the 2.0-10.0 keV unabsorbed flux 
of SWIFT J1749.4-2807 as measured from all available 
RXTE/PCA and Swift/XRT observations. Our dataset 
samples the last 7 days of the outburst, during which the 
flux decayed exponentially. We find that between MJD 
55306.5 and 55307.5 SWIFT J1749.4-2807 underwent a 
sudden drop in flux of more than an order of magnitude, 
less abrupt than the 3 orders of magnitude drop in flux 
obser ved in the previous o utburst of SWIFT J1749.4- 
2807 (jWiinands et alJ l2009). Similar dr ops in flux have 
been seen for other AM XPs (see, e.g. IWiinands et al.l 
120031 : iPatruno et al.l [2009a) . If we take into account 
the fact that SWIFT J1749.4- 2807 was first detected on 
MJD 55296 (|Pavan et al.ll2010l ). we estimate an outburst 
duration of about 12 days. 

8 http : //heasarc . gsf c . nasa . gov/Tools/w3pimms . html 
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2. — Pulse profile obtained by folding ; 
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Fig. 3.— Eclipse observed on observation 95085-09-02-02. The 
dataset starts on MJD 55302.9531, during the eclipse and shows 
that the egress occurs at orbital phase fb 0.282 (where orbital 
phase zero is the time of passage through the ascending node). 
We estimate the duration of the eclipse by assuming that eclipse 
is symmetric around orbital phase 0.25. The light curve was only 
corrected by the standard modeled background (see text). Count 
rates are in the 2-16 keV band. 

5. PULSATIONS 

Adopting a source position a — 17 h 49 m 31 s .94, 
6 = -28 o 08'05".8 (from XMM-Newton images, see 
IWiinands et al.ll2009D . we converted the photon arrival 
times to the Solar System barycenter (Barycentric Dy- 
namical Time) with the FTOOL faxbary, which uses 
the JPL DE-405 cphemeris along with the spacecraft 
ephemeris and fine clock correct ions to provide an ab- 
solute timing accuracy of 3.4 [is (jJahoda et al.l l2006h 

We created power spectra of segments of 512 sec of data 
and found strong signals a t frequencies o f « 517.92 Hz 
and « 1035.84 Hz (Alt amir ano et all 12010); these signals 
were not always detected simultaneously with a signifi- 
cance greater than 3er. 

To proceed furt her, we used the preliminary orbita l so- 
lution reported by IStrohmaver fc Markwardtl (|2010l ) and 
folded our dataset into 87 pulse profiles of « 500 sec 
each. We then fitted the profiles with a constant plus 4 
sinusoids representing the pulse frequency and its over- 
tones. We then phase-connected the pulse phases by fit- 
ting a constant pulse frequency plus a circular Keplerian 
orbital model. The procedure is described in detail in 
IPatruno et~all (|2010| ). In Table Q] we report the best fit 
solution and in Figure [5] we show one example of the 
pulse profile. 



It is known that the timing residuals represent a sig- 
nificant contribution to the X-ray timing noise, which if 
not properly taken into account can affect the determina- 
tion of the pulse frequency and the orbital solution (see, 
e.g. lHartman et~aTll2008t IPatruno et al]l2010T ). There 
is a hint of a correlation between the X-ray timing noise 
and the X-ray flux, especially between MJD 55302 and 
55303, where a slight increase of the X-ray flux is ac- 
companied by a jump in the pulse phas es, similarly to 
what w as reported for 6 other AMXPs by IPatruno et all 
(2009b). A complete discussion of timing noise in this 
source is beyond the scope of this paper and will be pre- 
sented elsewhere. 

In the middle panels of Figured] we show the fractional 
rms amplitude of the fundamental, and of the first, sec- 
ond and third overtone when the signal was > 3er signif- 
icant in ~ 500 sec datasets. The 95% confidence level 
upper limits are are estimated using ~ 3000 sec datasets 
(excluding detections) and plotted separately for clar- 
ity. When detected significantly, the rms amplitudes of 
the fundamental and 1st overtone are in the ~ 6 — 29% 
and ~ 6 — 23% ranges, respectively; the highest values 
are reached at the end of the outburst, where the un- 
certainties in our measurements also increase. Ampli- 
tudes for the fundamental as high as 15-20% rms have 
been seen before for at least one sou rce (although for a 
brief interval, see IPatruno et al. 2010), however, no other 
AMXP shows a 1st overtone as strong as we detect it in 
SWIFT J1749.4-2807. In order to compare the strength 
of both signals, in Figure Q] (lower panel) we show the 
ratio between the fractional rms amplitude of the 1st 
overtone and that of the fundamental. As can be seen, 
there are periods in which the ratio is approximately one, 
but also periods where the ratio is 2 or more. We note 
that these ratios are independent of the uncertainties on 
the background. 

6. ECLIPSES AND THE INCLINATION OF THE SYSTEM 

We searched the RXTE data for the occurrence of X- 
ray bu rsts and found none. Following iMarkwardt et all 
(2010), we also searched for possible signatures of eclipses 
and found two clear cases in the RXTE data (ObsIDs: 
95085-09-02-02 and 95085-09-02-04, beginning at MJD 
55302.97 and 55305.87, respectively). PCA data on MJD 
55306.97 (OBSID: 95089-09-02-11) samples an ingress, 
however, the count ra te is too low to extract useful 
information (but se e IMarkwardt fc Strohmaverl 120101 : 
iFerrigno et~al1l2011l ). 

The first and clearest case of an eclipse is shown in 
Figure [3] The average 3-16 keV count rate at the 
beginning of the observation is about ps 18.5 — 19.5 
cts/sec/PCU2 (only the standard modeled background 
has been subtracted) for the first ps 600 sec. Then the 
countrate increases within a few seconds to an average 
of ps 36 cts/sec/PCU2 and remains approximately con- 
stant for the rest of the dataset. The other dataset 
also shows a similar low-to-high count rate transition 
although at lower intensities: the observation samples 
less than 275 sec of the eclipse (at a rate of ps 18 — 19 
cts/sec/PCU2); the count rate after the egress is about 
ps 22 cts/sec/PCU2, i.e. much lower than in the previous 
case. 

Within the uncertainties on the unmodeled back- 
ground, both egresses occur between orbital phases of 



ps 0.2823 — 0.2825. During the eclipse the count rates 
in these two observations are consistent with the ex- 
pected background of ps 18 — 19 cts/sec/PCU2 , imply- 
ing that that SWIFT J1749.4-2807 most probably shows 
total eclipses; however, given the uncertainties in the 
background (see Section I3.1J) and the sensitivity of the 
PCA, this should be tested and better quantified with ob- 
servations from instru ments like XMM-Newton, Suzaku 
or Ch andra (see also iPavan et al.1 120101 : IFerrigno et al.1 

ImCD). 

Using the best solution reported in Table [1] we also 
searched for pulsations in the 600 second period during 
which the companion star is eclipsing the neutron star 
(see above). We found none. Upper limits are uncon- 
strained. 
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Fig. 4. — Inclination of the binary system vs. the the neutron 
star mass. For each point we also mark the mass of the companion 
star M c (in units of Mq) and the mass ratio q = M c /Mjvs- 



With our improved orbital solution and the measured 
times of the two egresses we determine the phase of egress 
to be no larger than 0.2825. Assuming the eclipses are 
centered around neutron star superior conjunction, the 
eclipse half- angle is ps 11.7°, corresponding to an eclipse 
duration of ps 2065 sec. 

We do not detect any evidence of absorption in the 
form of dips in the light curves, probably due to the 
fact that our dataset only samples ps 1.5 orbital peri- 
ods. These dips are common in other eclipsing LMXBs 
and thought to be due to the interaction of the photons 
from the central X-ray emitting region by structure on 
the disk rim or by what is left of the stream of incom- 
ing matter (from the companion) above and below the 
accretion disk. These dips are known to be highly en- 
ergy dependent; both eclipses and egresses in our data 
are energy independent. 

Assuming that the companion star is a sphere with a 
radius R equal to the mean Roche lobe radius, then the 
radius of the companion star can be approximated as 



0.49 • q 



2/3 



Rl ° ' 0.6 • g 2 /3 + i n (i + q i/3) ' W 

where a is the semi-major axis of the system and 
q = M c /Mns is the ratio betwe en the compani on and 
neutron star masses, respectively (|Eggletonlll983f) . From 



geometrical considerations in an eclipsing system, if the 
size of the X-ray emitting region is negligible compared 
with the radius of the companion star, then Rl is also 
related to the inclination i and the eclipse half-angle 4>: 



Rl 



■ 2 ■ 

sin i 



(2) 



when the eccentricity of the system is zero (see, e.g., 
Chak rabarty et al]|1993l ; also note that the half-angle 
of the eclipse might be smaller a s the star filling it s 
Roche lobe is not spherical, see, e.g.. lChanan et al]|1976| ). 
These two equations in combination with the mass func- 
tion form a system of equations that allow us to find 
the inclination of the binary system as a function of the 
neutron star and companion star mass. In Figure |4] we 
show our results. For a neutron star with mass in the 
0.8-2.2 (1.4-2.2) Mq range, we find inclinations in the 
74.4° - 77.3° (76.3° - 77.3°) range and companion mass 
in the 0.46-0.81 (0.67-0.81) M range. 

7. CONSTRAINING NEUTRON STAR PROPERTIES VIA 
PULSE-PROFILE MODELING 

Knowing the inclination to a high degree of precision 
is useful for pulse profile modeling to constrain neutron 
star properties including compactness and geometry. To 
explore what could be done, we tried fitting simple model 
lightcurves to the pu lse amplitude observa t ions (along 
the l ines explored bv iPechenick et all I1983L iNath et all 
120021 and lCadeau et al.ll2007D . The code we use has been 
tested ag ainst, and is in good agreement with, the results 
of lLamb et ail (|2009t ). 

We assume isotropic blackbody emission from one or 
two antipodal circular hot spots, and no emission from 
the rest of the star or the disc. At this stage we ig- 
nore both Com ptonization (which might be important 
Gicrlihsk i fc Pout anenT 2005) and disc obscuration. 

We consider as free parameters stellar mass and ra- 
dius, and the colatitude a and angular half-size S of the 



hotspot(s). Using only points where both fundamental 
and 1st overtone are detected with at least 3er signifi- 
cance, we search for models that fit all observations (am- 
plitude of fundamental and ratio of first overtone to fun- 
damental) and which have the same mass and radius. 
Hotspot size and position are permitted to vary between 
observations since accretion flow is expected to be vari- 
able. 

Although it is possible to obtain a high degree of har- 
monic content, due t o GR effects, from a single visi- 
ble hotspot (see also lLamb et~ai"1 l2009t ). we find that 
the strength of the harmonic is such that two antipo- 
dal hotspots must be visible in order to fit the data. We 
are also able to constrain system geometry. The la con- 
fidence contours restrict us to models with a ~ 50° and 
S = (45 — 50)°; the 2a contours permit a wider range of 
parameters but still require models where a = (40 — 50)° 
and 6 = (30 — 50)° (hotspots must be smaller if they 
are located closer to the pole). These results, within the 
frame of our simple model, suggest a substantial offset 
between rotational and magnetic pole in this source. 

Our models also put limits on stellar compactness. 
The la confidence contours exclude models with M/R > 
0.17MQ/km, while the 2a contours exclude models with 
M/R > O.18M /km. Although this does not rule 
out a ny common equation of state (Lattimcr & Prakash 
l2007f ). it does exclude some viable regions of dense mat- 
ter parameter space. 

Our simple calculations, while certainly not conclu- 
sive, illustrate the potential of this source. With better 
models, and phase-resolved spectroscopy using high spec- 
tral resolution observations, this system is an extremely 
promising candidate for obtaining tight constraints from 
pulse profile fitting. 
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